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Abstract: Four-dimensional colliding plane wave (CPW) solutions have played an impor- 
tant role in understanding the classical non-linearities of Einstein's equations. In this note, 
we investigate CPW solutions in 2n+2-dimensional Einstein gravity with a n+l-form flux. 
By using an isomorphism with the four-dimensional problem, we construct exact solutions 
analogous to the Szekeres vacuum solution in four dimensions. The higher-dimensional 
versions of the Khan-Penrose and Bell-Szekeres CPW solutions are studied perturbatively 
in the vicinity of the light-cone. We find that under small perturbations, a curvature sin- 
gularity is generically produced, leading to both space-like and time-like singularities. For 
n = 4, our results pertain to the collision of two ten-dimensional type IIB Blau - Figueroa 
o'Farrill - Hull - Papadopoulos plane waves. 
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1. Introduction 

Gravitational colliding plane wave (CPW) solutions have received much attention over the 
years, as a way to bring insight into the non-linearities of the collision of more realistic 
gravitational waves. The subject originated in the work of Khan and Penrose ||l| and has 
a vast literature, see for a review and an exhaustive list of references. 

CPW may play an important role in primordial cosmology as a possible seed for large 
scale structure formation or even setting the initial conditions in a pre-big-bang scenario 
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m. At a more formal level, colliding gravitational plane waves offer simple models in which 
to study the fate of the inner null singularity of realistic Kerr black holes [Q] , and a useful 
approximation to scattering at Planckian energies They also receive a cosmological 
interpretation as Gowdy universes |^] (i.e. with two commuting Killing vectors). 

While most of the CPW studies have taken place in the framework of four-dimensional 
Einstein gravity possibly with an Abelian gauge field, it is worthwhile to ask which of these 
results would continue to apply to higher dimensional gravity theories, and in particular 
to those which describe the low energy limit of string theory [^, ^, |9|, 11]. In particular, 



several maximally super symmetric plane wave solutions have been identified in type JIB |12| 
and 11-dimensional supergravity [|l^. It is an interesting problem to study their collisions, 
and whether spacelike singularities are generically produced in such processes. This may 
be especially tractable due to the high degree of symmetry of these backgrounds. Notice 
that in order to set-up the collision, it is customary to restrict to waves with bounded 
support along x"*", thereby breaking part of the supersymmetries. 

Finally, while exact CPW solutions of Einstein equations or their supergravity general- 
izations only pertain to very special initial conditions, it is important to study the stability 
of plane waves under small perturbations: indeed, plane waves usually come in an infinite 
dimensional moduli space, corresponding to the 2;^-dependent profiles of the various fields 
restricted by a single Einstein equation. A particle or string traveling along the x"*" direc- 
tion will generically be transmitted through the wave, but also be partly reflected back and 
alter the profile of the outgoing wave. This process is usually unaccessible in light-cone 
quantization approaches, as it involves the emission of = states. It is nevertheless of 
crucial importance when the background wave presents a null singularity, as it may turn it 
into a timelike or spacelike singularity. Studying this issue at the level of classical gravity 
may in particular shed light on the singularities observed in the closed-string parabolic 



orbifold [|14|, |l|] 



We start in section 2 by reviewing generalities on colliding plane waves in four dimen- 
sions, as well as several exact solutions which we will aim at generalizing. In section 3, we 
propose an ansatz (eq. ( |3.2D below) for 2n -|- 2-dimensional CPW solutions, which leads to 
same Einstein equations as in four dimensions, albeit different boundary conditions. Using 
this isomorphism, we construct some explicit solutions in 2n -|- 2 dimensions by upgrading 
four-dimensional solutions. In section 4, we develop a perturbative scheme which allows 
us to determine the solution in the vicinity of the light-cone for arbitrary boundary con- 
ditions. We use it to study the stability of gravitational and electromagnetic plane waves 
to small counter-propagating perturbations. Finally, we apply our method to construct 
higher-dimensional analogues of the Khan-Penrose and Bell-Szekeres metrics ^. We close 
in section 5 with a discussion of our results. The results of the perturbative computations 
can be found in Appendix A and B. 



^Another analogue of the Bell-Szekeres metric was constructed recently, but in the context of higher 
dimensional Einstein-Maxwell gravity jll]. 
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2. Colliding plane waves in four dimensions 
2.1 Generalities 

Colliding gravitational plane waves in four dimensions have been an field of intensive study. 
The metric ansatz compatible with the existence of two commuting spacelike Killing vectors 
dx,dy is given by the Rosen-Szekeres line element, 

cfs2 ^ 2e-'^Uudv + e-^ (e^ cosh Wdx^ - 2 sinh Wdxdy + e"^ cosh Wdy"^) , (2.1) 

where M, U, V, W are functions of the light-cone coordinates u, v only. We will restrict our 
attention to colliding waves with aligned polarization, ie = 0. In addition, one may 
allow for an electromagnetic field, 

F = i {dHi Adx + dH2 A dy) (2.2) 

where Hi and H2 are functions of u,v only. 

In studying plane wave scattering in flat Minkowski space, one usually assumes that 
space is flat ahead of each of the incoming plane fronts, say at u < and v < 0. Space-time 
is thus divided into four sectors: in the past region P : u < 0,v < 0, we have flat Minkowski 
space with M = U = V = W = Hi = 0; the right region R : u > 0,v < corresponds to 
the incoming left-moving plane wave, described by n-dependent profiles U{u),V{u), etc; 
similarly, the left region L : u < 0,v > corresponds to the right-moving plane wave, 
described by t;-dependent profiles U{v),V{v), etc. In the forward region F : u > 0,v > 0, 
the two waves start to interact, leading to a metric ( p.l| ) depending non-trivially on both of 
the light-cone coordinates. The problem is thus to determine the functions U, V, W, M, Hi 
in the forward region, given their values on the characteristics u = and v = 0. 

For this, notice that a change of u and v coordinates allows to set M = in the right 
(resp. left) region. The other functions U, V, W, Hi in the incoming region are freely chosen 
functions of u (resp. v), subject to the condition that the corresponding plane wave should 
satisfy the Einstein equation Ruu = (resp. Ry^ = 0)- Iii the interacting region, the 
Einstein equations require that is a free two-dimensional field: its value throughout 
region F is therefore determined immediately in terms of its boundary values at u = and 
v = 0: 

U = -log[f{u) + g{v)] (2.3) 

There is no loss of generality in assuming that /(O) = g{0) = 1/2. It is often useful then 
to change coordinates from {u,v) to {f,g)- Solving for the other functions V,W,Hi is in 
general a complicated non-linear problem, except for purely gravitational collinear waves 
{W = Hi = 0), where V satisfies a linear Euler-Darboux equation, and can be determined 
by Green's function techniques p^ , |2^ . When either W or Hi are non-zero, the problem 
is more difficult, although integr ability provides solution-generating techniques which in 
some cases allow to obtain exact solutions. Finally, M may be integrated by quadrature as 
it satisfies the free 2-dimensional Klein-Gordon equation with sources (Einstein equations 
for a higher-dimensional generalization of the ansatz ( |2.1[) will be displayed in Section 3, 

O-ii)). 
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Figure 1: Kinematical set-up for colliding plane waves. The dotted line denote a "fold" singularity 
in the incoming regions. 



While physically realistic waves have a smooth wave front, it is often useful to allow 
for mild singularities at w = or t; = 0, in order to describe idealized impulsive or shock 
wave profiles. As shown by O'Brien and Synge |jl^, the correct matching conditions to 
impose at n = are that the transverse metric gij and its derivative duQij be continuous. 
In terms of the ansatz (|2.1| ), this implies that M and V have to be continuous and f{u) 
and g{v) in ( |3.14| ) are at least across the boundaries u = and v = 0. In general then, 

f{u) - ^ ~ u^eiu), g{v) - ^ ~ -v^e{v), (2.4) 

where 0{x) is the Heaviside step- function and a, /? > 2. For a = 2, there are delta functions 



in components of the curvature tensors which have to be interpreted as distributions ||17|. 
Upon writing the incoming wave in Brinkmann coordinates, 

ds^ = 2dx+dx- + {H,^{x+)Xf + Hy{x+)Y^){dx+f + dXf + dY^, (2.5) 

it is easy to see that this corresponds to an impulsive plane wave, with H^^y oc 6{x^). The 
a = 4 case on the other hand corresponds to a shock wave, with H^^y oc 9{x~^). 



2.2 Exact four-dimensional CPW solutions 

Using the framework just outlined, many exact CPW solutions in four dimensions have 
been constructed over the years. We now briefly review several interesting solutions that 
we will be interested in generalizing to higher dimensions, see for an exhaustive review. 
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A two- parameter family of purely gravitational CPW solutions is given by the Szekeres 
solution [iril 



V 
W 
M 



—2ki arctanh 
H = 0, 



f 



+ 9 



2k2 arctanh 



9 



(2.6) 
(2.7) 



where / 



-2kik2 In 



u°' and g 




(21 



as in (pj) and A;f = 2(q - l)/a, k^ = 2{(5 - I)/ 13. A 



case of particular interest is the Khan-Penrose solution 
the incoming region i?, it corresponds to a profile 



which arises for a = (3 = 2. In 



2dudv + [1 + ue{u)fdx'^ + [1 - ue{u)fdy^ 



{2.{ 



which reduces to flat space both before {u < 0) and after (u > 0) the wave front. The 
singularity at n = describes an impulsive gravitational wave, with delta function profile. 
The singularity at n = —1 on the other hand is merely a coordinate singularity, often 
called "fold singularity" in the literature. The profile in the region L is identical to ( p.g| ) 
up to exchanging u and v. In the interacting region however, the geometry is curved, 
with a space- like curvature singularity at + = 1. The metric furthermore becomes 
complex at n > 1 or w > 1. It is thus legitimate to excise the region behind the space- like 
singularity in region F, as well as behind the fold singularity in regions L, R. Indeed all 
but exceptional u = cste causal geodesies coming from region P fall into the singularity in 

Another interesting explicit solution is the Bell-Szekeres solution [^], which describes 
colliding electromagnetic waves in Einstein-Maxwell gravity. The solution in the interaction 
region takes the very simple form given by 



ds'^ = 2dudv + cos^(u — v)dx^ + cos^(u -|- v)dy^ , 

A ' 



y 



sin(n + v) 



(2.10) 
(2.11) 



where the replacements u — > u6{u) and v v9[v) are implicit. In region R, the incoming 
wave may be written in Brinkmann coordinates, 



ds^ = 2dX+dX- + dX"^ + dY^ -\{X^ + Y^){dX 



(2.12) 



which we recognize the standard electromagnetic plane wave. One of the most interesting 
properties of the Bell-Szekeres solution is that the metric in the interaction region F is 



in fact diffeomorphic to a slice of the Bertotti-Robinson AdS2 x S'^ space [21, 22 1. In 
contrast of the purely gravitational collision described above, the collision of two collinear 
elecromagnetic waves is therefore free of any singularity. Instead, the Killing vectors dx 
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and dy become null at u + v = {2k + l)'/r/2 and u — v = {21 + 1)tt/2, respectively. There is 
therefore a Killing horizon at u + v = Tr/2, joining on to the fold singularities at u = 7r/2 
and V = 7r/2 in the incoming regions. In the literature, one usually excises the region 
behind the fold singularity and the Killing horizon, although this appears to be much less 
justified than in the Khan-Penrose- Szekeres case. 

Note that the incoming plane wave ( 2.12| ) is the four dimensional analogue of the 
ten dimensional type IIB plane wave found recently by Blau - Figueroa o'Farrill - Hull - 
Papadopoulos (BFHP) |]l2[. One of the goals of this note is to analyze the scattering of such 
PP-waves in higher dimensions, a situation which is of interest in string theory. In order to 
avoid undue suspense, let us immediately state that in contrast to the Bertotti-Robinson 
space, AdS^ x cannot be viewed as the collision of two BFHP plane waves. 

Finally, we would like to stress that there are two ingredients in the dynamics of 
colliding plane waves: firstly the equations of motions and secondly the matching of the 
solution in the interaction region F to incoming plane waves in region L and R and flat space 
in region P. Different matching prescriptions deflne different (un)physical problems. In 
particular, Gowdy-type cosmologies satisfy the same equations of motion (usually displayed 
using different coordinates) but different boundary conditions. 



3. Colliding plane waves in arbitrary dimensions 

Having reviewed the basic features and simplest CPW solutions of Einstein-Maxwell gravity 
in four dimensions, we now generalize these results to gravity in even dimension D = 2n + 2 
with a minimally coupled n + 1 form fleld strength, whose action is given by 

S = j A.V3^(fl--^F„..,„^.F« (3.1) 

Our main interest actually lies in the n = 4 case, which describes a subsector of ten- 
dimensional type IIB supergravity, upon restricting to self-dual conflgurations of the 5- 
form fleld strength. Other cases n = 2, 3 may also be relevant in the context of type IIB 
string theory and F-theory compactified on K^^ respectively. This motivates the following 
ansatze for the metric, 

ds^ = 2e-^Uudv + e'^^^'^Xdxl + ■■■ + dxl) + e'^^^+^Xdyl + ■■■ + dyl). (3.2) 

and for the n + 1-index anti-symmetric tensor fleld strength, 

FnA.1 = —= (duHidu A dxi A • • • A dxn + dyHidv A dxi A ■ ■ ■ A dxn 

+duH2du A dyi A ■ ■ ■ A dyn + d^H2dv Adyi A - ■ ■ A dy^^ . (3.3) 

Indeed, these ansatze are appropriate for incoming BFHP waves in type IIB supergravity 
||2^ (or their analogues in lower dimension), which preserve an SO{n) x SO{n) symmetry. 
Here M, U, V, Hi and H2 are all functions of u and v only. Factors of n have been inserted 
for later convenience. One may also consider non-collinear polarization, or switch on other 



-6- 



fields present in supergravity such as the dilaton, axion or other antisymmetric tensor fields, 
but we shall refrain from doing so. Our excuse is that, besides simplicity, the equations 
of motion satisfied by the ansatz ( |3.2| )-(^) are in fact identical to those satisfied by the 
four-dimensional Rosen-Szekeres solutions ( |2.1| )-( [2?^ )^. 

3.1 Equations of motion 

Indeed, just as in four dimensions ||T8| , |2|, the Einstein equations on the ansatz (|3.2| ) can 
be written as two chiral equations, 

dlu + djJduM - i- {{duUf + {duVf) 
dlu + d,ud,M - {{d^uf + {d,vf) 

and three integr ability conditions, 

dud^U - duUd-M = 0, (3.6) 

2dudyV - duVd^U - d^VduU + 2k {e^-^ duHid^Hi - e^+^ duH2d.,H2) = 0, (3.7) 

1 2n — 1 

dud^M - —duVd^v + — — djjd^u = o. (3.8) 

2n 2n 

Satisfying these equations automatically imply (|3.4| ) and (|3.5| ). In the L and R regions, the 
latter are simply the equations Ruu = and R^^ = satisfied by the incoming plane waves. 
Note that the constraint ( |3.4[) (resp. |3.5| ) implies that the mean radius of the transverse 
directions e~2^ is a concave function of u (resp. v) except when it vanishes. f{u) (resp. 
g{v)) is therefore a concave decreasing function, indicating that a (possibly coordinate) 
singularity in regions L amd R is inevitable. Thus "fold" singularities are generic in higher 
dimensions as well. Note also that these equations automatically imply the vanishing of 
the Ricci scalar 

R = -^e*^ [(2n + l)duUd^U + duVd^V - 2n(9„9„M + 2dudvV)] (3.9) 

Other curvature invariants i?2 = (Rfiu)'^ and i?4 = {Rf_iupa)'^ are however non-trivial, and 
as we shall see generically have singularities in the interacting region. Finally, it is easy to 
check that that the only conformally flat plane- wave solution satisfying the ansatz ( |3.2| )- 
( |3.3| ) is flat d + 2-dimensional Minkowski space, with the sole exception of the Bertotti- 
Robinson metric AdS2 x in 4 dimensions. In particular, AdS^ x does not fit into the 
ansatz (p.2|)-( |3^ ) in 10 dimensions. 

In addition to the Einstein equations (^.4|)~( |3^ ), the equations of motion for the n-form 
electromagnetic gauge field are given by 

2dudyH, - duHid^V - d^HAV = 0, (3.10) 
2dudyH2 + duH2d^V + d^H2duV = 0. (3.11) 

^This would remain true if we allowed non-parallel polarizations preserving the SO{n) x SO{n) symmetry. 



= ^(e^-V^(a,Fi)2 + e^+^(a„i72)^) (3.4) 
2n 

= ^{eU-y{d.H,f + e^^''{d.H2f) (3.5) 
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In dimension d = 6, 10 one can impose a self-duality of the fieldstrength F = *F which 
relates the two functions Hi and H2 by 

duH2 = e-^duHi, d,H2 = -e-^d,Hi. (3.12) 

This is notably the case in type JIB supergravity. In dimensions d = 4, 8 there exist no 
self-duality condition in Lorentzian signature, nevertheless ( |3.12| ) may still be imposed. For 
simplicity, we shall impose ( p. 12 ) for all values of n, and denote H = Hi. The constraints 
(|3j|)-(^) thus reduce to 

d'^u + dudM - — ({duf + (dvf) = -e^-^{dHf, (3.13) 



2n V J n 

where d stands either for du or 

As in four dimensions, ( |3.6| ) implies that e~^^'^''"^ is a free field in two dimensions, and 
therefore 

U = -log {f{u)+g{v)). (3.14) 

The functions f{u) and g{v) can be directly determined from the value of U on the char- 
acteristics n = and v = 0, assuming /(O) = g{0) = 1/2. For purely gravitational waves 
{H = 0), Eq. ( ^ ) is the linear Euler-Darboux equation, and can be solved by the same 
Green's function techniques that applied in the four-dimensional case. When H ^ 0, the 
equations ( |3.7| ) and ( p.lO| ) are non-linear, but identical to the four-dimensional case: they 



can thus be dealt with just the same integrability techniques. The constraints (3.4)- 
can then be integrated to yield the remaining function M. Changing coordinates from u, v 
to f{u),g{v), they reduce to 

a,M . l±±(a,Vf . y^S.Hf - '^^^ + . 0. (3.15, 

,,M + l±li,0,vy- + te-^-(0,Hf - + . 0. ,3.10) 

As in 4 dimensions, the inhomogeneous terms may be absorbed by defining 
2r? — 1 

M{u, v) = ^— log(/ + g)- logiduf) - log{d,g) + Sifiu),giv)). (3.17) 
The constraint equations then simplify to 



dfS + - (fc e-^ (dfHr + ^{dfVy ) = 0, (3.18) 



dgS + -[ke--'id,Hr + ^{d,Vr) =0, (3.19) 

involving derivatives with respect to /, g only. 

This shows that the equations of motion for the 2n -|- 2-dimensional ansatz (p.2|)-(p 
with "self-dual" n + 1 field-strength are isomorphic to those for the n = 1 four-dimensional 

v3 



collinear ansatz (|2.l| )- (l2.2| ), upon identifying' 



Un = Ui, Vn = Vi, Hn = Hi, (3.20) 

1 n — 1 

M„ = -Ml + ( log if + g)- log(9„/) - log{d,g)) . (3.21) 



^This isomorphism remains valid upon dropping the assumptions of collinearity and self-duahty. 
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This immediately allows us to construct 2n + 2-dimensional solutions from four dimensional 
ones, as we now discuss. 



3.2 Higher dimensional CPW from four-dimensions 



As discussed in the previous section the form of the solution has a very simple dependence 
on n. In particular in the "self-dual" case the equations for U and V are independent 
of n. This implies that a CPW wave solution in d = 4 (n = 1) automatically gives a 
solution in d = 2n + 2, with n > 1. A word of caution however is that, due to the non- 
homogeneous terms log{duf), log{dvg) in ( |3.21 ), the properties of the resulting solution 
can change significantly for n > 1 . In particular the existence of a smooth matching of the 
solution for n = 1 does not in general imply a smooth solution in higher dimensions. In the 
following we will analyze the properties of higher dimensional solution for some particular 
four dimensional solutions. 



Purely gravitational higher dimensional CPW 

To illustrate this construction, we start with the four-dimensional Szekeres solution (| 
with 

f{u) = -- u'^e{u), g{v) = -- v^e{v) (3.22 



Using the upgrading procedure outlined above, one finds a four-parameter family of purely 
gravitational CPW solutions in arbitrary dimension 2n + 2 



U 
V 

M 



- ln(l — u 
-2ki arctanh 



1 - vl^ 



2k2 arctanh 



1 - u° 



- In duf + ^ log(l/2 - /) - \nd,g + ^ log(l/2 - g) 
2n 2n 

^ 2n-l-(fci+fc2)^ _ „a _ + M in(l -u^) + ^ ln(l - v^) 
2n 2n 2n 



+ 



kik2 



n 



ln(l - - t;^ + 2u"v'^ + 2^^i"v/3(l - u")(l - z;/^)) . 



(3.23) 
(3.24) 



(3.25) 



where we implicitely replace u" — > u'^9{u) and — > v^6{v). As in four-dimensions how- 
ever, the two parameters A;i,/c2 get related to a,/? by demanding that M be continuous 
across n = and w = 0, 



a — 1 



kl = 2n , kl = 2n , 



a 



(3.26) 



The family of physical solutions relevant in dimension 2n + 2 is therefore different from 
that relevant in dimension 4. The fact that a curvature singularity arises at ua + = 1 
remains nevertheless true. The solution satisfies the appropriate junction conditions across 
the null surface u = (resp. u = 0) if a > 2 (resp. (3 > 2). In region i?, the incoming 
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plane wave metric is given by 

ds^ = (1 - u^eiu))'^ 2dudv + (1 - u'^e{u)) " „^ ^ 



1 — u2 0{u) 



+ (1-«°«W)- (3.27) 

In order to transform this metric in Brinkmann coordinates one has to redefine the u 
coordinate for ti > 

2n-a 

(1 - u"e{u)) = du. (3.28) 
Then the metric can be brought into the form 

ds^ = 2dudv + e^iufdxj + el{u)dyf (3.29) 
and the standard change of variables into Brinkmann form gives 

ds^ = 2dx+dx- + {H^{x+)Xf + Hy{x+)Y^){dx+f + dXf + dY^, (3.30) 

where 

1 d'^Cx 1 d'^ey 

An analogue of the four-dimcnsional Khan- Penrose solution may be obtained by setting 
a = P = 2, i.e. kf = fe^^: the resulting incoming profile 

1 71 1 /y»H~ 

= -j=6{x+) + ^7-— — (3-32) 

yn ^2 (1 — [x+Y) 1 

Hy = -^S{x+) - ^ '''^ e{x+). (3.33) 

ni (l-(x+)2)^ 

shows an impulsive (delta function) component, together with a non-vanishing tail which 

depends on x^: this is markedly different from 4-dimensional case, where space was flat 
on either side of the wave front. We shall return to the true analogue of the Khan- Penrose 
solution in Sections 3.3 and 4.2. 

Similarly, a higher dimensional analog of the Szekeres solution is obtained by ki = 



k2 = ^ 2n^^^i^ which gives a = f3 = m with integer m > 3 and one finds 

= hx{x+)e{x+), Hy = hy{x+)e{x+), (3.34) 

where the behavior of hx, hy as x"*" — > depends on the value of the integer m. For m = 3, 
hx,hy diverge as x"*" 0. For m = 4, hx,hy has a discontinuity at x~^ 0, but also 
exhibit a x^ dependent tail: this is in contrast to the 4-dimensional Szekeres solution, 
which corresponded to a true shock wave with constant Brinkmann parameter on either 
side of the wave front. For m > 4 the functions hx,hy vanish as x''" — 0. 

To summarize, by upgrading the general Szekeres solution we found exact purely grav- 
itational SO{n) X SO{n) symmetric collinear CPW in higher dimensions. They describe 
different incoming profiles from the 4-dimensional case, yet exhibit the same features of 
fold singularities in the i?, L region and space-like singularity in the forward F region. 
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Electromagnetic CPW solutions 

We now upgrade the four dimensional Bell-Szekeres solution ( |2.2| )-( prTl| ) into a higher 
dimensional colliding plane wave solution with flux. Solving (|3.4| ) and (|3.5| ) determines M 
to be 

n — 1, f cosiu — v) cos(u + v)\ 

M= log ■ - ^9 ^ ' 3-35 

n \ sm[2u) sm(2v) J 

H = = ^=sm(u + v) (3.36) 
V2k 



hence the metric is of the following form 

( n-l 
sin(2ti) sin(2v) 1/ \, , 2 \ 5./ \, , 2 \ 
7 r -. r- auav + COS" (n + ulfaxi +...) + cos" [u — v)[ay^ +...]. 
cos(m — V) cos(u + v) J 

(3.37) 

Unfortunately, this metric exhibit a curvature singularity on the characteristics n = and 
u = 0, as can be seen by computing the scalar curvature invariant 

2(n-l) 

cos(ti — v) cos(u + v) ' 



- i sin(2^)sin(2.) J ^^'^^^ 

This implies that the solution in the interaction region cannot be glued smoothly to in- 
coming plane waves in region L and i?, and is therefore unphysical. Another way to reach 
the same conclusion is to change variables to n = (sin u) " ,v = (sin v) « , so as to 
remove the singular part of the dudv term in ( |3.37| ): then / = 1/2 — ■u2"-i , g = 1/2 — -02"-! 
have exponents a, (3 smaller than two, in contradiction with the O'Brien Synge conditions. 
Finally, note that the incoming plane waves in the L and R regions are markedly different 
from the BFHP plane wave in higher dimensions, which was our motivation to look at the 
Bell-Szekeres solution in the first place. 

While this attempt to produce a physical higher-dimensional CPW solution from the 
Bell-Szekeres solution has failed, there are nevertheless four-dimensional electromagnetic 
CPW solutions which lead to acceptable higher-dimensional ones. For example, recall that 
using solution- generating techniques developped in four dimensions, a purely gravitational 
CPW solution Uq,Vq,Mq may be turned into an electromagnetic one by the following 
transformation : 

U = Uq, y = Vb - 2 In ( cos^ a + sin^ a g-^o+^o) (3 39) 

^sinacosa fe-^o+^o - 1) 



k cos^ a + sin^ a e-'^o+^o 



M = Mo-21n(cos^a-hsin2ae-^o+^«), H -- 

Starting from the four-dimensional purely gravitational Szekeres solution ( |2.8| ), one may 
thus obtain a 5-parameter family of higher dimensional CPW solutions with flux. A special 
choice of parameters then leads to a solution which is regular at the characteristics u = 
and f = 0, although it displays a curvature singularity at f{u) +g{v) = in the interaction 
region. 
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3.3 Toward higher-dimensional Khan-Penrose and Bell-Szekeres CPW 

The construction in the previous section, based on upgrading known four-dimensional so- 
lutions into higher dimensional ones, has failed to produce what deserves to be called 
higher-dimensional analogs of the Khan-Penrose and Bell-Szekeres solutions, namely solu- 
tions describing the collision of two purely gravitational impulsive plane waves, and two 
purely electromagnetic shock waves, respectively. This must simply mean that we failed 
to impose the appropriate boundary conditions on the characteristics. 

Indeed, a higher-dimensional analog of the Khan-Penrose solution should describe the 
collision of two plane waves with Rosen coordinate metric 

ds"^ = 2dudv + [1 + u9{u)f{dx'l + ■■■ + dxl) + [1 - ue{u)f{dyl + ■■■ + dyl) , (3.40) 

flat on either side of the wave front at n = 0. The appropriate boundary conditions are 
therefore 

U{u, v) = - log[(l - u^)" + (1 - v'^f - 1] (3.41) 
y (n, 0) = n log[(l + u)/{l- u)] (3.42) 
V{<d,v) = n\og[{l + v)/{l-v)] (3.43) 

A solution could be obtained using the Green function technique developed in 
however for n > 1 the resulting integrals appear to be too difficult. 

Similarly, a higher-dimensional analog of the Bell-Szekeres solution should describe the 
collision of two BFHP plane waves with Rosen coordinate metric and F field 

ds^ = 2dudv + cos^ u {dxj H h dxl + dyj A h dy^) (3.44) 

f'U 

H = 2 cos^" udu (3.45) 







The appropriate boundary conditions are therefore 

U{u, v) = - log[cos2" u + cos^" t; - 1] (3.46) 
V{u,0) = , V{0,v) = (3.47) 

H{u, 0) =2 cos^" udu , F(0, v) = 2 cos^" vdv (3.48) 
Jo Jo 

Again, it is not clear how to obtain such a solution with the known solution-generating 

techniques in four-dimensions. In the following we will follow a different approach, using 

perturbation theory around the light-cone. 

4. Perturbative plane wave collisions 

In this section, our aim is to analyze higher dimensional plane wave collisions in a pertur- 
bative expansion around the light-cone - or, equivalently, when one of the waves is of much 
smaller amplitude than the other. As an application, we shall obtain approximations to the 
higher-dimensional Khan-Penrose and Bell-Szekeres CPW solutions as defined in section 
3.3, in the vicinity of the light-cone. We will also be able to study generic perturbations 
of the impulsive and shock plane waves in higher dimensions. 
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4.1 General set-up 

We consider an incoming left-moving plane wave in 2n + 2 dimensions, given by the Rosen 
coordinate metric and flux 



ds'^\ji = 2dudv + e 



e^y''H-){dxl + ■■■ + dxl) + e-Ty''^-\dyl + ■■■ + dyl) 



Fn+l = ^duH^^\u) du A {dXi,..n + dyi...n) (4.1) 

Here we took advantage of coordinate reparametrization invariance to set M = 0, and 
restricted to self-dual flux configurations. The functions U^'^\V^^\ H^'^^ are assumed to 
satisfy the chiral equation of motion ( p. 41 ) . We also assume that these functions vanish at 
u < 0, and satisfy the appropriate regularity conditions at u = 0. 

On the other hand, the incoming right-moving plane wave is assumed to be of the form 



ds^lr = 2dndt; + e-^^'''(™) 



e^^'^'(-)(dx2 + . . . + dxl) + e-^^'^'(™)(rfy? + • • • + dyl) 



Fn+i = ^duH^^\ev) du A {dxi,..n + dyi,„n) (4.2) 

where we introduced a parameter e which can be viewed as the relative strength of the 
two plane waves: upon going to Brinkmann coordinates, the Brinkmann mass parameter is 
proportional to e. As for the right-moving wave, we assume that vanish at 

V < 0, and have at most an impulsive singularity at f = 0. For simplicity, we also restrict 
to collinear polarization. 

We shall be interested in a perturbative expansion in e, therefore in the vicinity of 
the characteristic axis v = in the forward region F. To this purpose, we expand the 
left-moving wave profile as 

5(1) > 0) = i + ^g2V^ + ^gsv^ + ... (4.3) 

> 0) = ViV + ^V2V^ + -^^V3V^ + ... (4.4) 

H^^) {v>0) = hiv + ^h2v'^ + ^hsv^ + ... (4.5) 

where fl'^^^ (?;) is defined as usual by C/*^^^ (w) = — log[5^"^^ (?;)-|-l/2]. The vanishing of (71 , , /lo 
follows from our assumptions on the singularity at f = 0. The chiral equation ( |3.5| ) allows 
to eliminate e.g. g^^^ in favor of V^^\ H^^^ order by order in v, e.g. 

g2 = -^{hj + vi^) , g3 = ^{2hih2-hivi + 2viV2) (4.6) 

54 = [hf (1 - 2n) + vi^ + 4/iin2 {h^ - 2h2Vi) (4.7) 

+hl ((n - 2) (2n - 1) vi^ - 2n^V2) + n {-3vi^ + 4n (/i2^ + + ^1^3))] 

Our goal is now to solve for the solution in the interacting region. Using the field equation 
( |3.4| ), U is determined throughout the interacting region by its value on the characteristics 
uv = 0: 

U{u, v) = - log f/(°) + g^^^ (ev)] (4.8) 
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where we defined as usual U^^\u) = —log[f^^\u) + 1/2]. The other functions V,H,M in 
the ansatz ( |3.2| )-( |3^ ) can be determined order by order in e by Taylor expanding 



H = fW(^x) + evHi(u) + ieV/72(n) + ^eVi^sN + ... 
V = + evViiu) + ^eVFaH + ^e^v^Vsiu) +... (4.9) 



M= + et;Mi(u) + ^e2t;2M2(u) + ie3v3M3(n) +... 



The validity of this expansion may be checked by self-consistency, or by looking at the 
known Khan-Penrose and Bell-Szekeres solution for n = 1. 

4.2 Higher-dimensional Kahn-Penrose solution 

We now specialize to the case where the incoming right-moving plane wave is a purely 
gravitational impulsive profile ( |3.4CI| ), hence 



U{u, v) = — log 



(4.10) 

At leading order in e, the equations ( |3.6| )~( |3^ ) for Hi{u),Vi{u), Mi{u) are first order 
homogeneous linear ODE's, 

H[ - = (4-11) 

1 — 

VI - -^V, = (4.12) 
1 — ti^ 

M[ - rr^Vi = (4.13) 
1 — ti^ 

Using the boundary conditions Hi{0) = hi, Vi{0) = vi, Mi(0) = 0, we obtain 

Hi{u) = hi{l + u)"/2 (1 _ uyn/2 (4_;^4) 

Vi{u) = VI {1 - u^y/^ (4.15) 
Mi(n) =^;in2FiQ,l + ^,|n2^ (4.16) 

where the hypergeometric function reduces to a simple algebraic function when n is odd. 
The norm of the Killing vector dx,y thus becomes 



\d,J^ = (1 - uf ± -vi{l - u)\l - u^)-'^/\ev) + 0{e^) (4.17) 
n 



which implies that the Killing horizons are shifted to 



\dx,y\'^ = : evr^^n 2V2a"/2 (4.18) 



where A = 1 — u. This is in fact an upper estimate, since higher order terms in ( 4.17] ), 



if singular at u = 1, may lead to a critical exponent higher than n/2. At this order, the 
curvature invariants R2 and R4 remain zero however. 
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At second order in e, the equations for H2{u),V2{u), M2{u) become first order linear 
ODE's with a source, 

H'2 - -^Hi = hivmil - n)-"-i (4.19) 
\ — 2, 

^2 - Y^^i = \H-{n + 1)/.?] (1 - u^)—' (4.20) 

^2 - Y3^^2 = -^^^(1 - tx2)-"-i [(2n - l)hl + {n- l)vl] (4.21) 

These equations are too complicated to give a general answer for any n, however explicit 
solutions can be found in Appendix B for the values n = 1 ... 4 of interest. The main 
result here is that the curvature invariant i?4 is now non-zero^ and in fact diverges in the 
interacting region at u = 1, as i?4 ~ A"+^. 

We have computed the solution up to order for n = 1 ... 4 in the special case where 
the left-moving {v) profile is purely gravitational, i.e all hi = 0. It is then consistent to 
set H[u, f ) = throughout the interaction region. This allows us to obtain the curvature 
invariant i?4 to second order in e (see Appendix A for explicit results). We find in general 
that higher order corrections diverge faster. Looking at the zeros of l/i?^, we find that 
higher order corrections shift the precise location of the pole. For example, in the n = 4 
case relevant for type IIB, we find 



Rl' = ^2 + -T^i^^r + 0{e') (4.22) 

3vf ovi 192 

which shows that the pole is shifted to^ 

i?4 = 00 : ev = 0(A"+2) (4.23) 



Determining the precise coefficient would require to resum the series in (4.22). In fact, more 
divergent terms in higher order contributions to could even increase the exponent, so 
that the relation provides only an upper estimate. To fourth order, we find that the critical 
exponent for the location of the Killing horizon ( [4.18| ) remains equal to n + 2. In contrast 
to the n = 1 case however, the numerical coefficients appear to be different: it is therefore 
conceivable that the space-like singularity may be hidden behind a Killing horizon, in 
agreement with cosmic censorship. 

Finally, a case of particular interest is when the two colliding waves have identical 
impulsive profile: this is the higher dimensional analog of the Khan-Penrose solution. The 
solution may be obtained to order by setting 

vi = 2n , V2 = , V3 = An , V4 = , (4.24) 

in the formulae in Appendix A. 



''Since they arise only at order e^, we systematically rescale R2 and R4 by a factor of e^. In this purely 
gravitational case the Ricci square R2 remains in fact zero to all orders in e. 
^This relation in fact holds for any n to order e''. 
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4.3 Higher-dimensional Bell-Szekeres solution 

We now consider a purely electromagnetic shock wave as our right-moving background 
profile given by ( 3.44 ). The function U is thus determined throughout the interaction 
region by 



U (n, v) = — log 



cos^"n + g^^^(ef) - ^ 



(4.25) 



For the most part, we shall consider purely gravitational left-moving perturbations, i.e. 
Vi = 0. We shall however relax this assumption in the n = 1 case when discussing the 
stability of its horizon. 

To leading order in e, the equations ( ^ ) ^(3^) then reduces to the linear system 

Vi = — ^ , h'; + n^Hi = (4.26) 
n cos" u 

Using the boundary conditions Vi(0) = 0,Hi{0) = hi, one therefore obtains 

Sill TIU 

Hi{u) = hi cos nu , Vi(u) = -hi (4.27) 

cos" u 

One may check that to this order, the curvature invariants R2 = {R^u)"^ and = {R^upa)'^ 
remain zero. The length of the Killing vectors dx,y however is shifted to 

\\d.,yf = e-(^^^)/" = cos^^ ± ^^^ehi + 0{e^) (4.28) 

n cos" ^ u 

For n = 1, the first order correction vanishes at u = vr/2, hence the horizon remains to this 
order at u = ti/2. For higher n however, the correction term blows up at u = ti/2. Were 
the second order term not more singular than the e order term, the horizon would be 
shifted to 

hiev ^ ~ even (4.29) 

hiev ~ n — — j ~ n odd (4.30) 

Unlike the Khan-Penrose case, it turns out that higher order term do change this behaviour, 
and lead to higher exponents. 

At order e^, the equations of motion become linear with a source, 

^ H'2 ^ 3cos[(2n-l)^]+cos[(2n + l)u] ^2 (4 3^) 

n cos" u 8 cos^"+^ u ^ 

n 2„ {n - 1) sin[2nn] + 2n sin[2(n - l)u] ^ . . 

H^+n H2 = -n — h^ (4.32) 

4 cos"+^ u 

, nsin[(2n — l)nl — (3n — 2) sinn , 9 

M!^ = ^ , 9 \ hi (4.33 

^ 4ncos2"+iu ^ ^ ^ 

These equations are readily integrated for given values of n, using the boundary conditions 
-^2(0) = ^125^2(0) = M2(0) = 0. The same pattern continues to hold to higher order, 
allowing us to get the solution to any order desired. Explicit solutions up to fourth order 
can be found in the appendix. We now briefiy summarize our results. 
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In four dimensions, the colHsion of two colhnear electromagnetic shock wave appears to 
be a rather smooth process. Indeed, when the two waves profiles are identical, the metric 
in the interaction region is given by the Bell-Szekeres solution, which is diffeomorphic to 
the Bertotti-Robinson or AdS2 x 5*2 metric. It does however present a Killing horizon for 
the two Killing vectors dx^y, at u + v = tt/2. For more general profiles, we find that the 
curvature invariants remain finite. 



R2 = Ri = Ahl + 8hih2€v + {hi + Ahl + 4:hih3)ievf + 0{e' 
To order e^, the norm of the Killing vectors is given by 

Il/l2/i2 



/j2 I 

\dxf = u^± hiuev + -jievf + -hih2{evf 



I92u 



To leading order, the horizon therefore lies at 

tt/2 = 



u 



±^hiev 



(4.34) 



(4.35) 



(4.36) 



As is well known however, the collision of an electromagnetic plane wave with a grav- 
itational wave does lead to a curvature singularity. Indeed, if we allow for a non-zero 
= viv + — „f ^ + . . . perturbation, the curvature invariants become 



R2 — Ra 



V 2 

.2 4(/ii cosu + sinti)((2/i2 — /ii^i) cosn + 2f2 sinn) os, 

4(/ii+t;i tann) k ev+0[e ) 

cos^ u 

(4.37) 

where we calculated but do not show the term. As u ^ 7r/2, the curvature blows up. 
The precise location of the singularity can be found by expanding I/R2 and keeping the 
dominant term at each order, 



Rr, 



4vl 



V2^ 

2vf 



-ev 



32A2 



lev 



(4.38) 



where we defined A = u — Tr/2. The curvature singularity therefore appears at v = ©(A^). 
For /ii = /i2 = this agrees with the explicit solution found by Griffiths p^]. 

In contrast, in 6 and higher dimensions, we observe that singularities are generically 
created, even for purely electromagnetic plane waves (y*-^^ =0)- To order e^, the curvature 
invariants have a pole of high order at u = it/2. The order contribution however is still 
more singular at n = '7r/2, implying that the singular locus is shifted away from u = tt/2. 
Similarly, the norm of the Killing vectors diverges at higher order in e, leading to a deviation 
from u = tt/2. Using the same technique as in ( 4.22| ), we may derive the upper estimates 
for the locus of the Ricci square singularity, Riemann square singularity and of the Killing 
horizon: 



R2 
R4 
\\dx,y\ 



00 



00 



: 



O 
O 
O 



u 



u 



TT^ 




2y 


n 


TT^• 




2. 




TT^ 




2/ 





(4.39) 
(4.40) 
(4.41) 
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Our results may be summarized in the following table of upper critical exponents for n = 1 
through n = 4, 



n 


a 


P 


7 


1 


n.a. 


n.a. 


1 


2 


2 


3 


2 


3 


4 


3 


3 


4 


4 


5 


4 



Generically, it therefore appears that, in contrast to the four-dimensional case, for 
n > 1 a line of curvature singularity is created tangentially to the light-cone at -y = at 
the fold singularities. One half of this line corresponds to a space-like singularity, while the 
other half is time-like. In order to determine the exact critical exponent of this singular 
line (instead of the upper bounds that we have derived), one would have to analyze the 
degree of divergence of the curvature invariants at u = 7r/2 to all orders in e. In particular, 
in the absence of such an analysis, an essential singularity at u = 7r/2 cannot be ruled 
out. It would be interesting to check whether these conclusions remain valid in the case of 
non-collinear {W ^ 0) or non-purely electromagnetic {V ^ 0) waves. 



Figure 2: Left: collision of two electromagnetic shock waves for for n > 1. Right: collision of two 
trains of opposite shock waves, for special choice of pulse duration. 

Finally, we may obtain a perturbative solution describing the collision of two BFHP 
shock waves in type IIB supergravity by specializing our analysis to n = 4 and H^^^ (v) = 
H(°\v), i.e. 

/ii = 2n , /i2 = , /i3 = -2n^ , /i4 = , . . . (4.43) 

The perturbative analysis performed above then applies in the vicinity of either of the light- 
cone axes, with the added simplification that all even coefficients h2,h4,hQ . . . vanish. It is 
therefore tempting to conjecture that the spacelike singularity emanating from u = 7r/2, v = 
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joins smoothly onto the one emanating from u = 0,v = tt/2, as depicted in Figure 
2. Furthermore, it is plausible that the time- like singularity emanating from u = {2p — 
l)7r/2 merges on to the space-like singularity emanating from u = {2p + 1)tt/2. Assuming 
this is correct, an interesting configuration arises when the incoming waves involves a 
succession of two shock waves of opposite amplitude separated by a critical time vr: using 
time reversal invariance, it is easy to see that the singular locus will form a closed line 
inside the interaction region. 

5. Discussion 

In this note we have studied colliding plane wave solutions of 2n + 2-dimensional Einstein 
gravity with a n + 1-form field strength. A natural ansatz ( |3.2| )-( p^ was shown to lead to 
the same equations as in ordinary four-dimensional Einstein-Maxwell gravity. By upgrad- 
ing known four-dimensional solutions, we were able to construct exact solutions in 2n + 2 
dimensions. In general however, the boundary conditions suitable in dimension 4 are not 
appropriate in higher dimensions. In the second part of the paper, we therefore devel- 
oped a perturbative scheme in order to determine CPW solutions with specified boundary 
conditions, in the vicinity of the light-cone. 

The general conclusion of this study is that, just as in four dimensions, space-like sin- 
gularities usually develop in the interaction region, emanating from the fold singularities 
in the incoming plane waves. An exception is the case of four-dimensional purely elec- 
tromagnetic plane waves, which only develop a Killing horizon in region F. Switching on 
an arbitrarily small gravitational perturbation however immediately leads to a space-like 
singularity, as was already observed in the context of Kerr black holes Q . In higher dimen- 
sions, purely electromagnetic plane wave collisions already display space-like singularities, 
as well as Killing horizons. By specializing to identical incoming waves, we were able to 
obtain the metric of the higher dimensional version of the Khan- Penrose and Bell-Szekeres 
CPW solutions in the vicinity of the light-cone. For n = 4, the latter describes the col- 
lision of two BFHP plane waves in type IIB string theory, or rather of their shock wave 
generalization. 

As in earlier studies of plane wave collisions, it is worth stressing that the creation of 
singularities for arbitrary small wave amplitude is tied to the plane wave symmetry, and 
would presumably acquire a threshold in the case of more realistic pp- waves. 

Another general result is that the space-like singularity generally appears in combina- 
tion with a time-like singularity, the two of them meeting at the fold singularity. In earlier 
studies of the four-dimensional problem, it has generally be assumed that the region behind 
the fold singularity should be excised, thereby getting rid of the time-like singularity in 
the forward region. However, in view of recent developments in the understanding of the 
asymptotic structure of plane waves and holography |^, ^] , this prescription seems 

hard to justify. Surely, in the absence of the counter-propagating wave, the fold singularity 
is simply an artefact of the Rosen coordinate system. Even in the presence of the second 
plane wave, there exist null geodesies which can cross the fold singularity and possibly 
reach the time-like singularity without encountering any horizon. This does not necessary 
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pose a conflict with cosmic censorship, as the latter may not apply in non-asymptotically 
flat geometries. 

In contrast to pp-wave backgrounds, colliding plane wave geometries in string theory 
receive a' and gg corrections. It would therefore be very interesting to find exact conformal 
field theory descriptions of such geometries. The Bertotti- Robinson AdS2 x 5*2 geometry 
or the Wess-Zumino-Witten model 5/(2) x SU{2)/M. x M are exact solutions of string 
theory j^, but it is not known how to impose the junction conditions at = while 
preserving conformal invariance. On the other hand, string probes propagating in a pp-wave 
background are just another example of plane wave collisions, and it is a very important 
problem to understand their backreaction on the plane wave background. 

In this respect, let us note that there exists a different perturbative scheme than the 
one considered in section 4, where one expands in the amplitude of the perturbation g^^^ {v) 
rather than its gradient, i.e. define U{u,v) = —log [/'^"^ — 1/2 + e'^{g^^\€v) — 1/2)] . At 
each order in e the profiles are general functions of (n, v) rather than of u only. Assuming 
that the perturbation has compact support along v, one may compute perturbatively its 
effect on the background right-moving wave. To leading order in e, the equations (|3.6D-(pI^) 
are total derivatives with respect to v; after integration, they become ordinary differential 
equations with respect to u only, with initial conditions set by the left-moving profile 
at every value of v. This reduction to a one-dimensional dynamical system is 
very reminiscent to the dimensional reduction that takes place near a space-like singularity, 
however the instability of null singularities render this observation less useful. 

To conclude, plane wave collisions are an inexhaustible source of space-like singularities. 
It would be very useful to develop a holographic description of them, as it may provide 
some insight into the dynamics of strings near a cosmological singularity. 
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A. Perturbative expansions — the impulsive gravitational case 



A.l n = 1 

At second order, 



H2{u) = h2\ + hivi — A.l 

V 1 — u 2(1 — u) 



At third order, for vanishing hi, 

Vs{u) = (A.4) 

2(1-^2)2 

u((l + u^) vi^ + SuVl — v?viV2 + 4 (1 - V?') ^3 ) 

M3(u) = ^ , ^ (A.5) 

4(1-m2)2 

The fourth order correction was computed but is not displayed here. It allows us to compute 
the curvature invariant i?4 through order e^, 

Gvi^ 3uVl - v?vi^ + 2v2 - 2u'^V2 

= n ^ + ^'""^ n (A-6) 

(1 — u'') (1 — ■u^j 

^ ^2 (-1 + 43^2) vf + 55txvT^^i;2z;2 + 6 (l - u"^) + 8 (l - u'^) V1V3 

+3(e^^) — ;-4 + C(e ) 

2(1 — 'u^) 

Its inverse can be expanded around u = 7r/2, keeping the most dominant term at each 
order in e. 

The space-like singularity therefore lies at 

ev = C»(AV2) (A_g^ 

Similarly, the length of the Killing vectors reads, keeping the most dominant term at each 
order in e. 



2 t;iA3/2 v2A3/\ ,2 i;fAi/2 



,4 



The Killing horizon thus lies at 
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A.2 n = 2 



For conciseness, for n > 2 we restrict to the purely gravitational case, hi = 0. At second 
order in e. 



V2{U) 



2 {AV2 + U {Vl^ - 4UV2)) - (1 - u^) vi^ (log i^) 



8(-l + n2)^ 



(A.ll) 



^'^""^ = 64 (-1 + .2) ^) 

We computed but do not display the third and fourth order contributions. This allows to 
compute the curvature invariant R4 through order e^. 



2{3 + u^)vl 



R4 



2\ ,.2 



+24 (-3 + 2u^ + n^) V2 + 9 (-3 + 2u^ + w^) vf log j 
Expanding the inverse around 7r/2 



(A.13) 



1 + u 



u 



R, 



.1 2A^ 5A^ . 87 



64 



we find that the curvature singularity lies at 

ev = 0(A) 

The length of the Killing vectors reads 

11 ^ 1,9 .9 ^1 A %V2 — f? , / n9 v^ , 



hence the Killing horizon lies at 



-MevY - 
64 ^ ' 384A 

ev = 0(A) 



6144A2 



(A.14) 

(A.15) 

(A.16) 
(A.17) 



A.3 n = 3 

At second order. 



V2{U) 



M2{u) 



u 



(-3 + 2v?) 



+ 



V2 



u {uvl (3 - 6«2 + 2u'^) + 12v2 (-3 + 2u'^) (1 - n2)3/2) 
36(l-n2f 



(A.18) 



(A.19) 



We computed but do not display the third and fourth order contributions. This allows to 
compute the curvature invariant R4 through order e^. 



i?4 



2{2v? + 3)^;^ 
(1 - n2)5 '^^'3(l-n2^7 



uVl^^ (89 - 36tx^) v\ (A.20) 



-18(l-u2)^(3 + 2u2)v2 +0(e2) 
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Expanding its inverse, 
we find a curvature singularity at 



ev = C>(A3/2) ^^22) 
The length of the Killing vectors reads 

\\d,f = A2 - ^^^^ev - ^{evf ^^^M' + ..oI'oL. M' (A.23) 

" 6V2 5184V2A5/2^ ^ 248832A4^ ^ ^ ^ 



hence the Killing horizon lies at 

ev 



A. 4 n = 4 

At second order, 

-2u (-5 + 3n2) + 32(1 - ^2)^2 + 3(1 - u'^)^vi'^ log(4±^) 
V.(») - > ^3,(,_i.). ^ ^ 

4m (-23 + 42u^ - 39tx4 + 12^^) vi^ - 32(l - u^f (-5 + 3^2) ^2) 

M2(u) = J ^ (A.26) 

^ ^ 1024(1 -u2)^ ^ ^ 



3 (4u (-5 + 3u^) v\ - 64(1 - u^)'vi + 3(l - u^)^^ {x^^ 1^^^ 1^ 

1024(1 - u^f 



We computed but do not display the third and fourth order contributions. This allows to 
compute the curvature invariant i?4 through order e^, 

= + [-2- (-53 - 18n^ + 27n^) (A.27) 

+ (1-^x2)2(1 + ^2) j^96^2 + 27Vlog(^^)) +0(e2) 
Expanding its inverse. 



1 16A6 llt;iA4 41w2a2 

- — ^ -I ^ 

3uf 6i;i 192 

we find that the singularity lies at 

tv = 0(A2) (A.29) 

The length of the Killing vectors reads 

\\Q^f = a2 - "^ev - -^{evf + ^^{evf + '^-^^—{evf (A.30) 

" " 16 512A^ ' 8192A3^ ' 2621440A5 ^ ' ^ ' 

hence the Killing horizon lies at 

ev = 0(A2) (A.31) 
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B. Perturbative expansions - the electromagnetic case 

B.l n = 1 

At second order, 

1 



H2{u) = /i2COsn— ^hlsmu , V2{u) = —h2tanu , M2{u) = (B.l) 



At third order, 



3 

Hz{u) = h^cosu — -hih2smu (B.2) 

Vsiu) = -J^^ [5hl + 4/i3 + {hi + 4/13) cos 2u\ (B.3) 
o cos^ u 

M^iu) = ^/ii/t2tan^u (B.4) 



At fourth order. 



H^iu) = I [(-5hlh2 + 8/14) cosM - (3hf + 12hl + 16/11/13) sinu + 5/if/i2sec«] (B.5) 
8 

VJu) = (45/t?/i2 + 8/14 + (7hlh2 + 8/14) cos 2u) (B.6) 

Idcos^ u 

M4(u) = l(hi + 2hl + 4/11/13) tan^ u (B.7) 
8 

The Ricci square reads 

i?2 = Ahl + Shih2ev + {h\ + 4/ii + 4/ii/i3)(eu)^ + ©(e^) (B.8) 
One may check that this is also equal to the Riemann tensor square i?4 = i?2- 

B.2 n = 2 

At second order, 

H2{u) = /12 COS 2n — /i^ tan M cos 2?x (B-9) 

V2{u) = -2/i2tanti (B.IO) 

M2M = /i2^^(5cos2u- 1) (B.ll) 

At third order, 

, , (9/i2 (cosu + cos3u) + 2/i? (3sinu — 2sin3'u)tan('u)) , 
Hsiu) = /^3 cos 2u - hi^^ ^ ^ ^ B.12 

DCOSM 

y3(u) = ^ (4 {hf - 3/13) - hlsec{uf (6 + sec^ ix)) tan-u (B.13) 

= 3/^1/^2(1 + 3 cos(2^))tan^n 
16 cos"' u 
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At fourth order, 

H4{u) = hi cos2n + ^^:^M [4 {^h\ - 9hl - Uhihs) cos(2 



+hf (-72 + 40sec(ii)^ + seciu)"^^ + /i?/i2 (-14 + 101 cos(2u)) t&n{u) 



M4{u) 



tanu 



^2 (58/i?/i2 - 24^4 - hih2 (43 + 36 cos(2ti)) sec'' u) 
sec^(ti) — 1 



1152 



(529/if - 1296/i^ - 1872/n/t3 



(B.15) 

(B.16) 
(B.17) 



+ {-815hi + 432/1^ + 432/11/13) sec(u)^ + /tf (179 + 17cos(2u)) sec(n)^) 
The Ricci square curvature invariant reads 



i?2 = 6 [{hi + h2evf + hih^{ev) 



cos^ 2u 



(B.18) 



COS 2?/. 

-^i(e^)^^:5 6~ (1 - 90 cos 2u - 31 cos 4u) tan^ u + ©(e^) 

COS li 



In order to find the location of the singularity, we expand I/-R2 around u = 7r/2, keeping 
the dominant term at each order in e 



(B.19) 



At this order, the Ricci curvature singularity therefore lies at v = C'(A^). As explained in 
the text, this is only an upper estimate, as higher order e-corrcctions to R2 may be equally 
or more singular at u = 7r/2. We now turn to the Riemann square curvature invariant. 



i?4 = 



(B.20) 



/if (7 + 3 cos(4'u)) sec(u) 
2 

+/11/12 (7 + 3 cos(4u)) sec(ti)^ev 

+^{vef \l2 (-31/1^ + 48/ii + 48/ii/i3) + 12 {bbh\ - 48/i^ - 48/11/13) sec(u)^ 
48 L 

-5 (139/if - 48/ii - 48/11/13) sec(u)'' + 604/ifsec^ u - 263/ifsec^ u + 90/if sec^° nl + 0{e^) 



At each order in e, the most dominant terms are 



(B.21) 



The Riemann square singularity is therefore at v = C>(A^). The length of the Killing vector 
reads 



||a.,,r = A^±/ilAe.-^(e.)^+^^^^3 

The Killing horizons are therefore at f = O(A^). 



(B.22) 
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B.3 ri = 3 

At second order, 



H2{u) = /i2 cosS-u — ^/tf sec^ ^(sinn — sinS-u + sin5n) (B.23) 
V,{u) = -h^'-^ (B.24) 

COS'^ u 

O 

M2(u) = hi— ^(-3 + 56cos2u + 43cos4u) (B.25) 

^oo COS U 



At third order, 



H^iu) = /i3 cos 3^ — /''^^^^^ r360/i2cos^ u (2 cos 2u — 1) 
240 cos*' u '- 

+h\ (-80 sin u + 95 sin 3u - 89 sin bu)] 
V'i{u) = tanu [464^f - 240/i3 + sec^ u (-836/i? + 6O/13) 

+/i?sec^ u (450 - 103sec^ u + lOsec^ u)] 

, , / N /ii/i2 (33 + 104 cos 2u + 67 cos 4-^) tan^ u ,^ , 

-''^al'") = ^ TT-. 7 (B.26) 

^ ^ 144cos4u ^ ' 

At fourth order, 

Hi(u) = h4 cos 3u + 216/ti^2 <^°^ + [(-4312/ii/i2 + 120^4) cos 3^ - UOh^ cos 4u 
-48 (536/if - 45/ii - GOhihs) smu+ 16 (l76/if - i5hl - 60/11/13) sin3'u 
-138/if sec^utan?x + sec^M(8359/i?/i2 - 10396/if tanu) 
+sec^ u (-807/1? /i2 + 1927/it tanw) 

-12 sec M (2430/1? /i2 + (-2156/if + 45/ii + GOhihs) tanu)] (B.27) 

V4{u) = — (l79/i?/i2 - 15/14) + TTTT^ — [-10072/i?/i2 + 120/14 (B.28) 

J.0 i^U COS u 



+hlh2 (5400sec^ u - 1181sec^ u + 95sec(u)^) 



Mi{u) = ^ ^g^QQ^^ [4 (31049/if - 15075/1^ - 20475/ii/i3) (B.29) 

+sec^ u (-255964/if + 36900/i^ + 47700/ii/i3) 
+sec(ii)^ (178196/if + 450/i^ - 900/11/13) 

+9/ifsec(u)^ (^-6301 + 1184sec(u)^ + 50sec(u)^) 

The Ricci square reads 

i?2 = 8 \{hi + h2evf] _ ^hiievf (3sec2« - 4) [-2252/1? + 72O/13 

'- cos* M 45 

+sec2 u (3849/if - 54O/13) + 2/i?sec^ u (-585 - 362sec^ u + Ulsec^ «)] 

+C»(e3) (B.30) 

Keeping the most singular terms at 7r/2, we find 

1 47 

R-^ = — -^A^ + \A evf + OU^) (B.31) 
2 72/i? 6480A2^ / V ; v ; 
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The Ricci singularity is therefore at v = C(A^). The Riemann square reads 



/j2 

i?4 = h— (8 + 2 cos 2u + 7 cos 4u + 2 cos 6u + cos 8u) (B.32) 

3 cos*^ n 

2h\h2 , , ^ ^ 

H 3 — ( — cos 2^ + 8 + 10 cos(4ii) + 2 cos 6« + cos 8m) eu 

3 cos* u 

^{evf [32 (-563/if + 180/ii + ISO/ii/ia) + 48 (923/it - 180/ti - 180/ii/i3) sec^ u 



135 

-6 (11261^^ - 1080/ii - I26O/11/13) sec^ + (92029/tf - 3510/ii - 4860/ii/i3) sec^ u 
-18 (4267/if - - 6O/11/13) sec^ + 30690/if sec^° u - 5040/if sec^^ ^ 
+420/ifsec^^'u] +0(e^) 



Expanding, 



i^ji = -^A^ - A4/i3A8 ei; - ^A2(e^;)2 + ©(e^) (B.33) 
4/i| (12 36 



It is singular at f ~ ©(A^). The length of the Killing vector reads 



Wd^yf = A2 ± h_^v + -^rle^)' ± -^M' + ^^M' + C?(e5) (B.34) 
II a;,aii 36A4^ ' 162A7^ ^ 3888A4^ ^ ^ ' ^ ' 

The Killing horizons are therefore at v ~ C>(A3). 



B.4 n = 4 

At second order, 



, , , , (sin -u — sin 3u) cos 4'u , 9 
iJ2(u) = /i2 cos 4u + 2^^ ^ h\ (B.35) 



COS"" n 
= , (sin.-sin3.) ^^ 

COS'^ M 

^ _ tan^ u(22 + 439 cos 2u + 442 cos Au + 249 cos 6u) ^ 

~ 2048cos6u ^ ^ ' ' 

At third order, 

Hz{u) = — g- (48/if + 10/i3 - 3 (24/i? - 5/13) cos 2n + 6 (8/1? + /13) cos (B.38) 

COS U 

—24hf cos 6u + hs cos 6u — \2h1h2 sin 2'u — 2Ah\h2 sin4n — \2h\h2 sin 6-^) 
y3(u) = 4 (-2 + sec^ u) {-Shi + /is + /i? (3 + 3 cos 2n + 2 cos 4n) sec^ u) tan n (B.39) 



M3(«) = j^^^^ g (278 + 823 cos 2u + 698 cos 4m + 377 cos Qu) tan^ u (B.40) 

J-U^4t COS 1i 
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At fourth order, 

COS ''\.U 

H4{u) = h^cosAu H 2 — (3/i2 (sinu — sinS'u)) (B.41) 

COS u 

COS ^i'iX 

+3/if ^ — (sinu — sinSn)^ + 2hihssec^ u (sinn — sinSu + sin5n — sm7u) 

cos^ u 

_^ /ti/t2tan u ^^^^ _^ ^2 cos 2u + 1038 cos 4u + 8 cos 6u + 505 cos 8u) 

56 cos^ u 

1/4 (ti) = 4/i4 tan n (tan^ u - l) (B.42) 

+ ^^in" (8 - 166 cos 2u + 64 cos 4n + 197 cos 6u + 8 cos 8u + 169 cos lOu) 
224 cos 10 ' 

,^ , , , , / 9 \ 505 — 519sec^'u + 121sec^'u — Tsec^u ,^ . , 

M4{u) = hihz (sec^ u-l) (B.43) 

,4 /196643 , 4 503857 o ^ 

+/if 384 sec^ u + 624 sec^ u - 528 sec^ u + sec^ u - 60 sec^° u 

^ V 2048 2048 

12 21sec^^w\ /i^ (278 + 823 cos 2-^ + 698 cos 4u + 377 cos 6u)tan2n 
+6sec^2u — — + 2^ ' 



512 ) 1024cos6tx 
The Ricci square reads 

i?2 = 10 [(/ii + h2evf + hih^{evf] (B.44) 

cos^ u 

cos^ 4ii 

-5/if (eu)^— — -TT- (-1141 + 1656 cos 2u - 2244 cos 4n + 1992 cos 6n + 761 cos 8u) 

4096cosi6„^ 

+0(e=^) 
Expanding 

It diverges at f ~ ©(A"^). The Riemann square reads 

77-96cos(2'u) + 48cos(4'u) + 5cos(8n),2 /t^ 
i<4 = 3 /ii (B.46) 

4cos^n 

101 - 144cos(2m) + 72cos(4tt) + 5cos(8tt) , , 

H 5 h,\ h'2 

2 cos" u 

+—{evf [320 (-761/if + 64/i| + Uhih^) + 640 (l273/tf - 64/ii - 64/^1/13) sec^ u 
128 

-16 [104105/it - 2752/i| - 3136/11/13) sec^?^ 
+16 {mi2lh\ - 2048/ii - 2624/11/13) sec^u 

- (3232777/if - 10304/t^ - I3376/11/13) sec^ « + 2103680/ifsec^° u - 714416/ifsec^^ u 
+108848/if sec^'' u - 5083/tfsec^^ u + 252/tfsec^^ ^i] 



Expanding 



r,-i 2 ,0 644/i2 63 , ,9 

^4 = TT^ev ^— r(ew)2 + 0(e^) (B.47) 

^ 113/if 1132/i3 8-1132A2^ ; V ; V ; 
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we see that the R4 curvature invariant is singular at v ^ C(A^). The length of the Killing 
vector reads 



The Killing horizons are therefore at ~ 0(A 



4\ 
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